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^§ ■ 1 Introduction 

in 

In paper pQ from the data on CP asymmetries in 5^(5^) — > p ± vr T , p + p 
^ ' decays and BABAR data on CP asymmetries in Bd(Bd) — > 7r + 7r~ decays we 

determine the value of angle a of the unitarity triangle: 

jg; a = 96° ±3° , (1) 

where only a tree quark decay amplitude b — > uudib — > uud) was taken 
^ ■ into account. The numerical values of angle a obtained from the considered 

decays are consistent with each other and with the value which follows from 
the global CKM fit. This observation testifies to the validity of a proposed 
approach. 

As the next step in the present paper we will study what changes in the 
values of a are induced by QCD penguins. Our aim is twofold. First, in 
this way we will get an estimate of the theoretical uncertainty of the value 
of a determined in paper 0. Secondly, we will get the formulas for CP 
violating parameters describing these decays which vanish when penguins 
are neglected (C p7T , A^p, C pp , C wn ). 

The angle shifts Aa we are interested in were estimated in paper [2]; 
however, in that paper FSI phases were neglected, that is why C pw = A P q P = 
Cpp = Ctt-k = follows from 2\ (the nonzero penguin amplitudes are a 
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necessary, but not a sufficient condition for C pn ... to be nonzero). We will 
take these phases into account. The asymmetries depend on the differences 
of FSI phases in the processes described by the tree and penguin diagrams. 
One source of these differences is an imaginary part of the quark penguin 
diagram, the so-called BSS mechanism of the strong phases generation [3] 
(see also paper 0J). The phase of the penguin diagram depends on the gluon 
q 2 which is transferred to uu pair, each quark of which goes to different tc ± - 
or p ± -mesons. In this way the value of q 2 depends on the light meson wave 
functions and we can estimate it only roughly. Another source of FSI phases 
is hadron rescattering and even less is known about the values of the phase 
shifts between the penguin and tree diagrams generated in this way. In view 
of this we will determine FSI phases from the experimental data on CPV 
asymmetries, and investigate to what values of a it will lead. 

In Appendix we present the weak interaction Hamiltonian which is re- 
sponsible for b — > uud transition and calculate the necessary matrix elements. 
Using these formulas in sections 2, 3, and 4 we study B — > pn, pp and im 
decays correspondingly and extract the values of angle a from the experi- 
mental data on CP asymmetries in these decays. We conclude in section 5 
with the averaged value of a and a general discussion. 

2 a from B d {B d ) -> ^tt ± 

The time dependence of the decay probabilities is given by jHj : 

dN(MB,^P^) = (1± ^ p)c -</ T|1 _ 9(CW±ACW)x 

x cos(Amt) + q(S p7r ± AS p7T ) sin(Amt)] , (2) 

where q — — 1 describes the case when at t = Bd was produced, while 
q = 1 corresponds to Bd production at t — 0. In the case of T(4S') — > BdBd 
decay the flavor of the beauty meson which will decay to pn is tagged by the 
charge of a lepton in the other beauty meson semileptonic decay. A partner 
decay starts clocks as well, r is Bd(Bd) life time, while Am is the difference 
of masses of (Bd, Bd) system eigenstates (it equals the frequency of Bd — Bd 
oscillations). 

From Eqs. ()AI5j) and ()AI7j) we obtain: 

M-+ = AV ub V* d [l-0.07e i{s - a) } = AV ub V* d [l -0.07sin5 + i0.07cos5] , (3) 
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M-+ = BV: b V ud , (4) 
a M — h A 

\- + = - — — = e 2ia -[l-0.07sin5±i0.07cos<5] , (5) 

P M + £5 

where parameters g and p enter the expressions for (B d , B d ) eigenstates and 
we have substituted a = tt/2 in the (small) second term in square brackets 
in Eq. (|A15jl . Analogously we get: 

M + - = BV ub V* d , (6) 

M+- = AV* b V ud [l - 0.07e i(a+S) ] = AV* b V ud [l + 0.07 sin 6 - 0.07z cos 5} , (7) 

x+ -_qM^_ =e 2i*B [l + 0.07 sin <5-i0.07 cos 5]- 1 . (8) 

From the expressions for the quantities C pn and AC pn j^]: 

± AC pw = j~|^' 2 (9) 



we obtain: 



AC, 



a 2 - b 2 a 2 1 + AC p7r 



~ a i ■ b 2 " b- 1 - AC, 



p7T 



fa/6) 2 

^ = 0.28sinJ (1 , (10) 



where a = \A\, b = |.B|. 

The Belle and BABAR averaged result for AC p7r is jU]: 

AC p7r = 0.22 ±0.10 , (11) 

which leads to: 

a\ 2 

= 1.56 ±0.33 . (12) 



b 

From the averaged experimental result [S] 

C pw = 0.31 ±0.10 (13) 

and Eq. (fTUj) we get: 

sin 5 = 4.6 ± 1.5 . (14) 
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We see that poor accuracy in the measurement of C pn does not allow to get 
any definite information on the value of phase S. 

The next observable we wish to discuss is CP asymmetry A P q P : 



\M+~\ 


2 - \M~+\ 


2 + 


M+~\ 


2 _ 


|M"+i 
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\M+~ 


| 2 + \M~+\ 


2 + 


M+~ 


P + 


M-+| 


2 



{a/bf 
1 + (a/bf 



0.14 sin 5 - \ ;/ , (15) 



which should be compared with the experimental result jSJ: 

Ag P = -0.102 ±0.045 . (16) 
From (HSJ) and (HHJ) we get: 

sin<5 = -1.2 ±0.5 , (17) 

and it differs from given in Eq. (jl4j) by 3.5 standard deviations. This is 
the largest discrepancy we encounter in this paper. Averaging these two 
numbers we obtain: 

sin 5 = -0.62 ±0.47 . (18) 

Finally we come to the discussion of the observables which are sensitive 
to the angle a: 

S P „ ± AS pn = ^Z\2 > ( 19 ) 

2a lb ~ ~ a 2 lb 2 - I 

Sp-K = T777T [sin 2a cos 5 — 0.07 cos S cos 5 — 0.07 - sin 5 sin 5} , 

1 + a 2 /b 2 a 2 /b 2 + 1 

(20) 

AS pn = ttt^I - cos 2« sin 5 ± 0.07 cos 5 sin 5 sin 2a — 

1 + a^/o^ 

a 2 /6 2 - 1 

- 0.07^/t- sin S cos S sin 2a] , (21) 

a 2 /b 2 + 1 

where the definition of the phase 5 is A/B = (a/b)e iS and in the small 
terms proportional to 0.07 in the expression for S pjT we have substituted 
cos 2a = — 1. 
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Let us start the analysis of the experimental data from AS p7T . According 
to 0: 

AS pn = 0.09 ±0.13 , (22) 

which is much less than one. According to Eq. (|12j) the factor which multi- 
plies square brackets in Eq. (j^UJ) (and in Eq. (}2*Tj) ) is very close to one, that 
is why it is the expression in square brackets which should be much less than 
one. The second and the third terms of this expression are really very small 
and we can neglect them. What concerns the first term, it is small when S 
is close to zero or tt: 

sin 5 = ASpn = 0.09 ± 0.13 , (23) 

where a small deviation of cos 2a from - 1 is neglected. 
Now everything is ready and from Eq. (J20j) we get: 

sin 2a = 5*^/ cos 5 + 0.07 cos 5 , (24) 

where we omit the last term in square brackets since it is negligibly small. 
In order to find angle a from the experimental data jHj: 

S pn = -0.13 ±0.13 , (25) 

we should determine the values of cos 5 and cos S. 

To propagate errors from sin S to cos S, we consider gaussian distribution 
for sin 5 truncated to physical region | sin 5 \ < 1 , transform it into (non- 
gaussian) distribution for cos 5 and take an interval containing 68% of prob- 
ability. 

In this way from Eq. (fTSj) we obtain: 

| cos<5| = 0.88 ±0.12 . (26) 

The average value of | cos<5| appears to be close to one due to the so-called 
Jacobian pick. 

Concerning 5 it follows from Eq. (}2*3*j) that |cos5| = 1 with very good 
accuracy. Depending on the values of phases 5 and 5 we get the following 
four domains for the angle a: 

5^0 , 5^0: sin 2a = -0.13 ± 0.13 + 0.06 

a = 92° ± 4° (27) 
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5^0, 5^n: sin2ct = -0.13 ±0.13 -0.06 

a = 96° ± 4° (28) 



5^tt, 5^0: sin2a = 0.13 ±0.13 + 0.06 

a = 84° ± 4° (29) 

5^7t, 5^tt: sin2a = 0.13 ±0.13 -0.06 

a = 88° ±4° . (30) 

Thus QCD penguins split values of a obtained without taking them into 
account: 94° -> 92°, 96°; 86° -> 84°, 88°. If BSS mechanism is valid, then 
only the domains given by Eqs.(|27j) and (J2T?|) remain (see also [7]). 



3 a from p+p 



The time dependence of CP violating asymmetry is described by the following 
formula: 

dN(B d ^p+pJ) _ dN(B d ^p+p~) 

acp(t) = rlNin dt + rr — dt + - = -C pp cos(Amt) + S pp sm(Amt) . 

dN(B d ^plp L ) dN(B d ->plp L ) rr 
dt dt 

(31) 

Let us remind that the longitudinal polarization fraction fi = 0.98 ± 
0.01 ± 0.02 [8J; f L = 0.95 ± 0.03 ± 0.03 and its closeness to one greatly 
simplify the extraction of CPV parameters from Bd(Bd) — > p + p~ decay data. 
These parameters are given by the following expressions: 

r _ 1 - l-VI 2 r- _ 2Im ^PP Zo 9 \ 

U W 1 i I \ 12 ' °PP i i I \ 12 ' 



1 ± A 



PP\ 



From Eq. (|A14|I we obtain: 

n M 1—0 (Mp 1 ^- 01 ) 

\ = l J l£E = e 2 - 1 ,^ , = e 2 "*(l + 0.14zsinae i5 ) , (33) 

pp P M PP 1 - O.OTe 1 ^) V ' ' V ' 

where we use the same letter 5 for FSI phase difference of the amplitudes 
generated by penguin and tree diagrams as in the case of B — > pn decays. 
These differences would be really the same if BSS mechanism dominates. 
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Comparing the averaged experimental result for C pp 

C pp = -0.03 ±0.18 ±0.09 0, (34) 

C pp = 0.0 ±0.30 ±0.10 0; (35) 

C™ p = -0.02 ± 0.17 (36) 
with the theoretical expression which follows from Eqs. (J32)) . (|33|) 

C pp = 0.14sinasin5 « 0.14sin<5 (37) 

we get: 

sin5 = -0.15 ± 1.2 ; | cos<5| = 0.88 ± 0.12 , (38) 

where the same procedure of error propagation as in the case of B —>■ pn was 
used and we sum statistical and systematic errors of sin 5 as independent. 
From Eq. (J32j) we obtain: 

Spp = sin 2a + 0.14 sin a cos 2a cos S = sin 2a — 0.14 cos 5 . (39) 

According to the recent measurements: 

S pp = -0.33 ±0.26 0, (40) 

Spp = 0.09 ± 0.43 0, (41) 

S e p f = -0.21 ± 0.22 , (42) 

and we get two domains for a: 

5^0: a = 92° ± 7° (43) 

5 w 7T : a = 100° ± 7° . (44) 

Just as in the case of B — > p7r decays only the first domain remains if |d>| < 
tt/2 0. 

Let us note that using the isospin analysis (which allows to prove the 
smallness of the penguin contribution) it was obtained: 

a = 100° ± 13° 0, (45) 

a = 87° ±17° 0, (46) 

where the error is mainly due to the uncertainty of the penguin contribution. 
Extracting this uncertainty and averaging last two numbers we get: 

a = 96° ± 7°(exp) ± ll°(penguin) . (47) 
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4 a from Bd(Bd) —> ty + ty 

The time dependence of CP violating asymmetry is given by the formula 
analogous to Eq. (J3T|) with the evident substitution of n instead of p. Eq. 
(J32*Jl with the same substitution is valid as well, while for the quantity X n7T 
from Eq. (jA13|) we obtain: 

= = e 2ia (l + 0.28z sin ae iS ) , (48) 

and what concerns letter S we should repeat the comment made after Eq. 
©. 

The experimental data of BABAR and Belle for CPV parameters S nn and 
C nn were controversial though at present (with the latest Belle results) the 
divergence diminishes. In view of this we will perform a two step analysis, 
taking at the beginning only BABAR results and then the averaged results 
of two collaborations. 

Comparing the theoretical expression 

C n7T = 0.28 sin a sin 5 w 0.28 sin 5 (49) 
with BABAR result 

C^ ABAR = -0.09 ± 0.15 (50) 

we get: 

sin 5 = -0.32 ±0.54 , (51) 

while for S n7r we have: 

S n7T = sin 2a -0.28 cos 5 , S^ BAR = -0.30 ± 0.17 . (52) 
From Eq. (}51~]) we get: 

| cos<5| = 0.9 ± 0.1 (53) 

and two domains of a corresponding to two signs of cos^ 1 : 

5^0: a = 91° ± 5°(exp) ± l°(theor) (54) 
Assuming \5\ < ir/2 |2] we would get only the first domain. 
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5 « vr : a = 107° ± 5°(exp) ± l°(theor) . (55) 

Belle result: 

C™ lc = -0.56 ±0.13 (56) 

deviates by 2.5a from BABAR and (if correct) would require considerably 
larger P/T ratio than we use in our paper. 
Finally, averaging (|5T)J) and (|5T)|) one gets 

C« = -0.37 ±0.10 , (57) 

and comparing with the theoretical expression ()49j) we obtain: 

sin 5 = -1.32 ±0.35 , (58) 

which leads to: 

|cos5| = 0.55 ±0.25 . (59) 
From the averaged experimental result |B]: 

SZ = -0.50±0.12(exp) (60) 

we get two domains 2 : 

5^0: a = 100° ± 4°(exp) ± 2°(theor) (61) 

5 « vr : a = 110° ± 4°(exp) ± 2°(theor) . (62) 



5 Conclusions 

We have analyzed CPV asymmetries in 5^(1?^) — > p ± ir T , P^Pl an( i 7r+7r_ 
decays induced by the charmless strangeless 6-quark decay b — > uurf. This 
decay can proceed through a tree or penguin diagram. As it was noted in 
PP when the penguin diagram is neglected, one obtains the values of the 
unitarity triangle angle a from CPV asymmetries in these decays which are 
consistent with each other as well as with the value of a which follows from 
the global CKM fit. However, in order to determine the theoretical accuracy 

2 Assuming |<5| < ir/2 |Jj which follows from BSS mechanism of S generation, we would 
get the first domain only. 
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of a extracted from the decays under study one should take the penguin am- 
plitude into account. This was done in the present paper, where the moduli 
of penguin over tree ratios were calculated with the help of the factorization: 

<M 1 M 2 \j l j 2 \B>=<M 1 \j 1 \BxM 2 \j 2 \U> , (63) 

while FSI phase shifts between the tree and penguin amplitudes were ex- 
tracted from experimental data. 

In order to determine numerical value of a one should average the values 
which follow from the considered decays. Since the phase shifts 5 can be 
different in B — > pp, tvtt and pir decays, we get too many possibilities. That 
is why let us limit ourselves to the theoretically motivated case \S\ < it/2. 

Averaging Eqs. and (|BT|) we obtain: 

5-0: a PPtmr = 98° ± 4° . (64) 
Averaging it with (j2*Tj) and (J2~9*j) we get two possibilities: 

abound = 95° ± 3° , or 

a b ^ md = 91° ±3°, (65) 

and the last one corresponds to the smallest possible value of a. In the case 
\5i\ > n/2, 5 = averaging Eqs.(|2%|). (JUJ), (p^j) we get the largest possible 
value: abound = 102° ± 3°. 

The global fit results for a are: 

«Ht = 94° ± 8° , «SL ter = 94 ± 10° . (66) 

Thus the accuracy of the present day knowledge of a can be close to that 
of/3: 

P = 23° ±2° . (67) 

I am grateful to A.V. Fedotov for great help in the treatment of the 
experimental data. I would like to thank A.E. Bondar and M.B. Voloshin for 
useful remarks. This work was partially supported by grants NSh-2328.2003.2 
and RFBR 05-02-17203. 
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Appendix 



The strong interaction renormalization of the tree Hamiltonian which 
describes the beauty hadrons weak decays is much smaller than for the case of 
the strange particle decays since the masses of beauty hadrons are much closer 
to Mw in the logarithmic scale. In the leading logarithmic approximation 
for operators 0\ and O2 we have: 



G F 

Hi,2 = ~^VubV* d 



{ 


a 3 (m b ) 


4/6 




a s (M w )_ 





«7 a (l + 7 5 )k/7 a /l + 7 5 )w- 
-2/6 



+ 



a s {m b ) 



M M w) 

x d>y a (l + 7 5 )m + dj a (l + 7 5 )6w7 Q (l + 7 5 )w 



[m7 q (1 + 75)6 x 

(Al) 



and substituting as(Mw) = 0.12, as(mj) = 0.2, b = 23/3 we get: 

G 

H li2 = -J=V uh V* d {l.lu la {l + 7 5 )6d7„(l + 7 5 )m - 
- 0.2J7 q (1 + 7 5 )6m7 q (1 + 7 5 )u} . 



(A2) 



NLO calculations confirm and refine this result. From Table 1 of ^2] for 
the value A 4 = 280 MeV (which corresponds to as(Mz) = 0.118) we get 1.14 
instead of our 1.1 and -0.31 instead of our -0.2. 

At one loop the following QCD penguin operator is generated: 



H- 



G 



3-6 



^(v cb v; d + v ub v: d) ^^ 



, { M w 



In 



(d7 M (l + 75 )A6) x 



G 2 
x (u^Xu + dy„\d) = +_£y t6 y^0.03{--d7 a (l + 7s)& x 

x (uy a u + dy a d) + 2(d a 7 a (l + i 5 )b c )(u c -f a u a + d c y a d a )} , (A3) 

where A are eight colour SU(3) Gell-Mann matrices, and Fierz identity X a b^cd - 
-2/35 ab 5 cd + 25 ad 5 bc as well as unitarity relation V cb V* d + V ub V* d = -V tb V td 
were used. 

Substituting qy a q = \qy a (l + 7 5 )g + ^qj a {^ — 75)9 we find the renormal- 
ization factors +0.01 and -0.03 for operators 3 , 5 and O4, 6 respectively. 
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Cl = 


1.14, 


c 2 = 


-0.31 , 


c 3 = 


0.016 , 


c 4 = 


-0.036 


c 5 = 


0.010 , 


c 6 = 


-0.045 



At NLO the renormalization factors for the operators in which only the left- 
handed quarks are involved (O3, O4) are different from those for the operators 
in which both left- and right-handed quarks participate (05,Oq). From the 
same Table 1 of [TJ we get 0.016 and 0.010 instead of 0.01 and -0.036 and 
-0.045 instead of -0.03. 

Finally, the effective Hamiltonian which describes the charmless strange- 
less B d decays looks like: 

H = ^ [K«ciO! + c 2 2 ) - V tb V t * d (C3O3 + c 4 4 + c 5 5 + c 6 6 )] , 

(A4) 

O x = uj a (l + 75)6^70,(1 + 7 5 )w 

2 = dj a (l + 7 5 )&M7 a (l + 7 5 )u 

3 = dj a (l + 7 5 )6[m7 q (1 + 7 5 )u + dj a (l + j 5 )d) 

4 = d a i a {l + 75 )6 C [M c 7a(l + 7sK_+ 47a(l + Ts)^] 

5 = d7a(l + 75)&[m7«(1 - 75)" + dj a (l - j 5 )d] 
O e = d a j a (l + 7 5 )6 c [M c 7 a (l - 7 5 )u a + d c j a (l - 7 5 )d a ] 

(A5) 

and the complex conjugate Hamiltonian describes B d decays. 

Our next task is to calculate the matrix elements of H between B d and 
p ± 7r T , p + p~ and 7r + 7r~ states, which is the most difficult part of the job. 
We will present the matrix elements of 4-fermion operators as the product 
of matrix elements of two 2-fermion operators between Bd and a light meson 
and vacuum and another light meson. The validity of this factorization is 
questionable; in particular, in this approach the FSI phases due to the light 
meson rescattering vanish identically. However, we found the statement in 
the literature that the corrections to the factorization formulas are small, 
being proportional to A/m b or powers of as(rrib) [I3J 3 . In any case nowadays 
factorization is the only way to get expressions for the decay amplitudes from 
the fundamental Hamiltonian. 

Since we are interested in B d decays to charged mesons and we will fac- 
torize 4-fermion operators, let us present Eqs. (|A4j) . (jA5|) in the following 
form 

H = -^V ub V ud {axUj a (l + 75)^7^(1 + 7 5 )w - 



3 Sincc q 2 = (Ps d — P-k, p ) 2 = 0{m 2 r , mz) one can argue that the large distance contri- 
butions invalidate the factorization formula. 
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V ubV* d 

- 2a 6 u(l +7s)Kl 



a 4 M7«(l + 75)^7a(l + l5)U - 
7sH) , 



c 5 



(A6) 



-0.042 



where a% = c\ + |c2 = 1.04, 04 = 04 + |c3 = —0.031, a6 = C6 + | 
and Fierz identities ■07 a (l + 7 5 )y?X7<*( 1 + 75)^ = V>7«(1 + 7s)?7X7a(l + 75M 
^7„(1 + 75 VX7«(1 - 75)77 = -2^(1 - 75)w(l + 75 V were used. 

The matrix elements we are interested in were calculated in paper [2] 
assuming factorization. Up to a common factor which includes constant /„. 
and B — > 71 transition formfactor / ( m ^) f° r the amplitude of B d 
decay it was obtained: 

2m 2 n 

CL4 -\- ; — rd6 



7T + 7T 



M(B, 



TT ' 7T 



d 

VubV* d 



ai 



VubV* d 



(A7) 



(m u + m d )(m b 

where we use the result of [2] and take into account the difference of the 
rescattering phases of the tree (~ a±) and the penguin (~ 04 and ~ as) 
amplitudes 5. One evident source of this phase is the imaginary part of the 
penguin diagrams with intermediate u- and c-quarks. Let us demonstrate 
that only the last one should be taken into account in (|A7jl : 

M(B d -> n + n-)~V ub V: d (T + P(m u )) + V cb V; d P(m c ) + V tb V td P(m t ) = 
= V ub V: d [T + P(m u )-P(m c )}-V tb V td [P(m c )-P(m t )} . (A8) 



As we are interested in CP asymmetries we should calculate A^ 
TT+n-)/M(B d -> 7T+7T-): 



-2i 



?M(B d 



A 



-2i/3-2i7_ 



1 + 



PK)-P(mc) 
T 



VtbV t * d P(m c )-P(mt) 
VubV* T 



1 + 



P(m u )-P(m c ) _ V t * b Vtd P(m c )-P(m t ) 



V* h V ud 



l + (e 



2(0:— 7r) i(7r— a) 



sin 7 P(m c ) — P{m t ) 



(A9) 



sin /3 T 

where a, /3 and 7 are the angles of the unitarity triangle. 

Thus the absorptive part of P{m c ) contributes to 8. 

Since the penguin operator P(m c ) equals the correlator of two vector 
currents, one immediately picks up its imaginary part from the textbooks on 
QED. It depends on the gluon momentum transfer and when the square of 
this momentum transfer is much larger than 4m;?, we have: 



P(m c ) 



— ITT = — 



mi 



5 w 30° 



(A10) 
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Since u- and w-quarks to which gluon decays go to different light mesons, the 
value of the momentum transfer squared is determined by these mesons wave 
functions. It varies between and zero. Thus we see that the mechanism 
suggested in [3] leads to the small positive value of 5: 

S<30° . (All) 

If BSS mechanism determines the value of 5, it would confirm the validity 
of our approach. If, on the contrary, the large distance rescattering of light 
hadrons changes 5 substantially, one should await large corrections to the 
dispersive part (the ratio (P/T)) as well. In the present paper we will allow 
5 to vary between zero and 2ir, but we use the expressions analogous to (jA7|) 
for the decay amplitudes. 

Let us return to Eq. (|A7|) . With the help of the following equation: 

V: d V ub sin/3 y } 

and using the numerical values m u + m d = 11 MeV, m b = 4.5 GeV we obtain: 

M(B d - 7T+7T-) ~ V ub V* d l- e ^- a ^(-0.06)e^ 

sin p 

= V ub V: d [l + V.lAe^- a+5 ^ , (A13) 

where (3 = 23° and 7 = 63° were substituted (we are using the value of 7 from 
the global CKM fit in order to estimate a small correction to the amplitude) 
and we put a\ equal to one. In Section 4 we analyze the experimental data 
on CP asymmetries in B d (B d ) — > n + 7i~ decays. 

Coming to B d (B d ) — > p + p~ decays we should calculate the corresponding 
matrix element of the Hamiltonian presented in (A6). Factorizing 4-quark 
operators we observe that the term proportional to a 6 vanishes, since 
< p\d(l — 75)^)0 >= 0: the (pseudo) scalar current cannot produce a vector 
meson from vacuum. That is why instead of Eq. (jA7|) we get (see also |2]): 

M(B d -> pjp L ) V tb V td lS 

~ ai a^e , 

V u bV* d VubV* d 

M(B d — > PlPl) ~ V ub V: d {l + 0.07e^- a+ ^} . (A14) 

The production of transversely polarized p- mesons by the vector current is 
suppressed as (m p /ms) 2 , and the experimental data confirm the dominance 
of p L . 
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If BSS mechanism is responsible for the phase 5, then it should be the 
same as in Eq. (|A13|) . Eq. (|A14|) is used in Section 3 to extract the value of 
angle a. 

Our last problem is the calculation of the amplitudes M(B d — > p T n ± ) = 
M T± . The amplitude M h corresponds to p~ production from vacuum by 
(du) current, so the term proportional to does not contribute to it: 



M~ 



V ubV* d 



A 



a i 



VtbVt * d a a c lS - 

\r t/* 4 
V ubV ud 



AV ub V* H \l + 0.07e^ 7r - a+<5 - ) 



(A15) 



where A is the complex number. 

In the case of the amplitude M + ~ it is n~ which is produced from vacuum 
by (du) current, so the term proportional to a§ contributes as well: 



M+- 
VubV* d 



B{ ai - 



v ubV* d 



2ml 



a>4 



(m b + m u ){m u + m d 



-a 6 



(A16) 



see P]. Here B is the complex number. Unlike the case of B d — > 7r + 7r _ decay 
the terms proportional to and have opposite signs and as a result the 
expression in square brackets with good accuracy equals zero, leading to: 



BVaV. 



ub v ud , 



(A17) 



and the penguin pollution is absent (a\ is omitted since it is very close to 
one). 

The amplitudes of B d meson decays, M~ + and M +_ , equal to M + ~ and 
M~ + correspondingly with the complex conjugate CKM matrix elements. 
We will use formulas (jA15j) and (|A17|) in order to determine angle a in 
Section 2 and will omit index "— " from <5_, since 5 + did not enter Eq. (|A17|) . 

Thus the penguin pollution is minimal in pir mode, intermediate in pp 
mode and maximal in tttc mode (as it was noted in [2])- 
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